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We study the counterpart to the multi-photon down conversion in the quantised motion of a trapped 
atom. The Lamb-Dicke approximation leads to a divergence of the mean motional excitation in a 
finite interaction time for fc-quantum down conversions with k > 3, analogous to the situation in 
the parametric approximation of nonlinear optics. We show that, in contrast to the Lamb-Dicke 
approximation, the correct treatment of the overlap of the atomic center-of-mass wave function and 
the driving laser waves leads to a proper dynamics without any divergence problem. That is, the 
wavy nature of both matter and light is an important physical property which cannot be neglected 
for describing the motional dynamics of a trapped atom, even for small Lamb-Dicke parameters. 
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I. INTRODUCTION 

When the susceptibility of a medium interacting with 
an electromagnetic field of optical frequency depends 
strongly on the electric-field amplitude, one enters the 
domain of nonlinear optics. Nonlinear couplings of elec- 
tric fields of different frequencies usually emerge from an 
expansion of the susceptibility in terms of the electric- 
field amplitude. Prominent examples of such nonlinear 
couplings are second-harmonic generation or two-photon 
down conversion, which are due to a second-order sus- 
ceptibility X2- Nonlinear crystals have been successfully 
used to produce squeezed quantum states of light via 
a two-photon down conversion. The extension of two- 
photon down conversion to an arbitrary fc-photon pro- 
cess, where k > 2, has also been studied. While this 
might be viewed as a natural generalisation of the second- 
harmonic generation or the two-photon down conversion, 
it has been shown that there is a subtle problem in the 
theoretical description of such processes. Fisher, Nieto, 
and Sandberg ||] have argued that it is not possible to 
define states by applying the unitary time-evolution op- 
erator on the vacuum field state. This argument has 
been partially removed by a consideration using Pade 
approximants 0. Later on, however, it was shown by 
Elyutin and Klyshko || and Hillery @] that for k = 3 
and k = 4, respectively, a divergence occurs in the mean 
photon number for finite interaction times. This diver- 
gence property may be interpreted as an unphysical arte- 
fact coming from the improper treatment of the A-photon 
process. In fact, it has been shown that the usual para- 
metric approximation is incorrect in that it neglects the 
energy transfer and entanglement, between the pump and 
signal mode of the electromagnetic field which emerges 
when the pump mode is quantised [||-^). We note that 
the possibility of observing k = 3 nonlinear quantum opti- 



cal conversion processes in a damped cavity has recently 
been discussed ||. 

While these phenomena are well known and elaborated 
in the context of nonlinear optics, recent advances in 
laser-cooling pqjrj| , state preparation [^2 13 1, and detec- 
tion [jlij of the motional quantum state of single trapped 
ions, a new type of realisation of such nonlinear mode 
couplings became possible. Here the modes are repre- 
sented by the 3D harmonic center-of-mass oscillations of 
a single ion in the trap. The nonlinear mode coupling 
may be realized by appropriate laser irradiation which in- 
duces vibrational Raman transitions ||l5|-[l9[| . This opens 
possibilities to study such mode couplings with an almost 
perfect system (i.e. the motion of the trapped ion) where 
the damping of the motion is negligibly small apart from 
a small heating rate due to technical imperfections p0[ . 
In the Lamb-Dicke regime, where the atomic center-of- 
mass position is well localised with respect to the wave- 
lengths of the applied laser fields, mode couplings result 
which are analogous to the optical mode couplings in 
the parametric approximation. That is, a treatment of 
the dynamics based on the Lamb-Dicke approximation 
would reveal an unphysical divergence of the mean num- 
ber of vibrational quanta for fc-quantum processes with 
fc>3. 

Whereas for a trapped atom in the Lamb-Dicke regime 
one gets a close connection to the paramctrically ap- 
proximated optical couplings, for a trapped atom which 
is not well localised with respect to the laser wave- 
lengths, nonlinear modifications of the couplings occur 
which arise from the overlap of the atomic center-of-mass 
wavefunction with the laser waves, describing the mo- 
mentum transfer onto the atomic center-of-mass during 
laser- photon absorption and emission fl^Jjjfl . These re- 
coil effects may strongly influence the coupling strengths 
which depend on the number of excited quanta in the vi- 
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brational modes. They have been predicted j2l| and ob- 
served E| for the first time in the context of a nonlinear 
Jaynes-Cummings model, describing the dynamics of a 
laser-driven vibronic transition in the resolved-sideband 
regime. 

We will show in this paper that the nonlinear effects, 
caused by the overlap of light and matter waves, will re- 
move the unphysical divergence problem which arises in 
the Lamb-Dicke approximation. The paper is structured 
as follows: In Sec. II we introduce the effective Hamilto- 
nian for the motional dynamics of the trapped atom and 
we briefly discuss the justification and the validity of the 
parametric approximation for optical couplings and the 
Lamb-Dicke approximation for the motional couplings. 
The laser-driven motional dynamics is then considered in 
Sec. Ill and the divergence problem in the Lamb-Dicke 
approximation is studied in Sec. IV. In Sec. V the more 
realistic treatment of an unspecified degree of localisation 
of the trapped atom is shown to remove the divergence 
and some examples for the time evolution are given. A 
summary and some conclusions are found in Sec. VI. 

II. MOTIONAL COUNTERPART OF 
MULTI-PHOTON DOWN CONVERSION 

For the fc-quantum vibrational-mode coupling, we con- 
sider here a two-photon vibrational Raman transition 
which has been experimentally realised Jl^Jl^] and theo- 
retically studied in the context of nonlinear couplings of 
vibrational modes |l5|-[l^]. By application of two laser 
beams which are off-resonant with respect to a strong 
electronic dipole transition and which are detuned rel- 
ative to each other by multiples of the vibrational fre- 
quencies in the trap, vibrational Raman transitions can 
be driven which may be used to realize a quantum me- 
chanical counterpart of nonlinear optics, see Fig. [y. For 
an appropriate laser-beam propagation geometry which 
affects only the dynamics in one vibrational mode of fre- 
quency v, in the rotating-wave approximation the Hamil- 
tonian describing the effect of the Raman laser drive on 
the dynamics of the vibrational mode is given by jl5| 

H k = hnf k (a^a;ri) (ir/a) k + h.c. (1) 

where a and fit arc the annihilation and creation oper- 
ators of vibrational quanta, respectively. Here the laser 
difference-frequency has been chosen to be kv, i.e. fc 
times the vibrational frequency of the mode with k > 
1 p2| |. The effective two-photon coupling strength is 
given by k, and 

" =2 *—x — - 

is the so-called Lamb-Dicke parameter describing the lo- 
calisation of the spatial extension of the center-of-mass 
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FIG. 1. fc-quantum motional coupling by application of 
two off-resonant laser fields with laser difference-frequency 
= kv, where v is the frequency of the vibrational mode 
which is specified by the beam-directions, ki— k2, of the beat 
node of the two lasers. 

wavefunction in the ground state relative to the wave- 
length A of the beat-node of the two laser beams. The 
Hermitian operator functions fk{a^a;ri) strongly depend 
on rj and are defined by the normally ordered expressions 

A(ata;,) = e-W 2 ^(_iy_^L_^ (2) 

= : (2ryVata) * J k (r)V¥aj 2/2 : . 

These nonlinear operator functions correspond to over- 
lap integrals of the motional states \n), \n + k) and the 
beat-node of the laser fields. They account for the re- 
coil effects during the process of absorption and emis- 
sion of laser photons of the trapped atom. Since they 
depend only on the quantum number h = a* a, in the 
basis of its eigenstates, n\n) — n\n) (n — 0,1,2,...), 
these operators are diagonal, with their diagonal ele- 
ments fk(n; rj) = {n\f k {a*a] rj)\n) being given by 

^(-^)=(^)i4 fe) (^)e-" 2/2 , (3) 

where L^\x) are the associated Laguerre polynomials. 
For a well-localised atom, that is, for very small Lamb- 
Dicke parameters 77 <gC 1 , or more precisely for a small spa- 
tial extension of the atomic wavefunction r\\Jn + 1 <C 1 , 
one reaches the so-called Lamb-Dicke limit. Here usu- 
ally the Lamb-Dicke approximation is made, which takes 
into account only the lowest-order terms in r;. In our de- 
scription of the fc-quantum coupling ([l]) the Lamb-Dicke 
approximation is performed by replacing the operator- 
valued function fk(aJa; rj) by its limiting value for a small 
Lamb-Dicke parameter, 

lim/ fe (n;?7) = (4) 
77— >0 fc! 
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By replacing in the Hamiltonian (Q) the operator func- 
tion f k (n; rf) by the c- number given in Eq. (|]) one obtains 
the interaction Hamiltonian of the fc-quantum pro- 

cess in the Lamb-Dicke approximation, that is, in lowest 
order of the Lamb-Dicke parameter, 



HI ' = TiK k a" 



hut 



it* 



(5) 



with 



(irf) 1 



being the /c-quantum coupling strength in the Lamb- 
Dicke approximation. 

In the context of nonlinear optics the Hamiltonian (|^) 
describes the A-photon down conversion process where a 
is the signal mode and the pump mode has been paramet- 
rically approximated by replacing its operators by ampli- 
tudes whose values are included in the coupling strength 
Kk. It is well known that for k > 2 the parametric ap- 
proximation described by the Hamiltonian (|^) leads to a 
divergent behaviour of the mean quantum number in the 
signal mode for finite interaction times, i.e. 



lim 

-ti— >At, 



(n(t 2 )) - (n(ti)) = oo, 



(6) 



for a defined interaction time 0< Aioo <oo. In nonlinear 
optics the parametric approximation of the pump mode 
fails due to pump depletion and the entanglement of sig- 
nal and pump modes which is essential in this type of 
interaction, leading to a energy conservation of the total 
number of photons in the pump and signal mode |l[] . The 
parametric approximation essentially neglects the energy 
transfer from the signal to the pump mode, leading to an 
unbounded increase of the energy in the signal mode. A 
quantum description of the pump mode is therefore re- 
quired, regardless of how strong the pump field actually 
is compared with the signal mode ||-|8| . 

For the case of a single trapped atom, the coupling- 
strength n k contains the classically approximated field 
amplitudes of the two Raman lasers. While in nonlinear 
optics the parametric approximation for the pump mode 
fails it is expected that for a Raman-driven trapped atom 
it is the Lamb-Dicke approximation, rather than the re- 
placement of the laser-field operators by their classical 
amplitudes, which leads to a divergent dynamics. Note 
that the validity of the Lamb-Dicke approximation is, 
in principle, in contradiction with a divergent motional 
excitation in the trap. Eigenstates of the trap potential 
should only be populated for r\\J n + 1 <C 1, that is for 
higher excitations the approximation is no longer valid, 
and a diverging mean excitation violates this require- 
ment. 

Therefore in all cases, even for a trapped atom with 
small Lamb-Dicke parameter (77 <C 1), we have to con- 
sider the full problem including the nonlinear operator 
functions f k {fi]rf) in the Hamiltonian (Q). As already 



noted, these operator functions describe the effects of mo- 
mentum transfer onto the atomic center-of-mass motion 
during the laser-atom interaction [ ]D||l9]j2jJ . They are of 
particular importance for higher vibrational excitations 
r\\Jn + 1 > 1 and they are discarded in the Lamb-Dicke 
approximation. 



III. EQUATIONS OF MOTION 

To study the time evolution of the mean quantum num- 
ber (njt)), we will start by deriving from the Hamilto- 
nian (|lj) the equations of motion for the populations of 
the vibrational levels. The Schrodinger equation is given 
by 



\m) = H k \m), 



where \ip(t)} is the state vector in the interaction picture. 
Using as basis the number states \n), we get the following 
equations of motion for the coefficients ip n (t) = (n\ijj(t)) 



dt 



-i[9k{n\rj)ij) n+k + gl{n-k;n) ip n „ k ] , (8) 



where g k (n; rf) is given by 



9k(n;v) = <iv) k \l^r^fk(n;v), (9) 

and g k (n; rf) = for n < 0. For notational simplicity we 
will omit here and in the following the time argument of 
ip n (t) and will only write ip n . 

The time evolution of the populations of the number 
states, P n = tp^tpn, are obtained from Eq. (j^) and its 
complex conjugate, 



dPn 

dt 



21m [g k (n;r]) 



g k {n-k;rj) C-fcV'n 



(10) 

To calculate the second time-derivative of Eq. one 
requires the time-derivatives of combinations of the type 
ipn4>n+k which are given by 



d_ 

Of 



ipn^n+k = igt(n; rf) (P n+k - P n ) (11) 

+ i [g k (n-k;7])ip^_ k tf; n+k - g k (n+k;n) ^*ip n+2k ] . 

Inserting Eq. ([ll]) into the time-derivated Eq. ([h]) the 
second time-derivative of the number statistics results as 



d 2 P n 
dt 2 



= 2\g k (n;n)\ 2 (P n+k -P n ) 
-2\g k (n-k;r 1 )\ 2 {P n -P n - k ) 



(12) 



2 Re 



g k {n; rf) g k (n + k; rf) ip* +2k ipn 
+g k (n-k\rf) g k {n-2k\rf) ip*ip n - 2k 
-2 g k (n; rf) g k (n-k; rf) V>* +fe V>n-fc 
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We are interested here in the temporal evolution of the 
mean quantum number 



d 2 (h(t)) ^ d 2 p n (t) 



dt 2 



n=0 



dt 2 



(13) 



which can be calculated with the help of Eq. (IT4). Here 
only the first two terms of Eq. ( p"2| ) contribute to the 
sum in Eq. (|l3|), whereas the real-part given in Eq. Jl2] ) 
cancels. The resulting equation of motion for the number 
statistics reads as 



d 2 (n(t)) 
dt 2 



OO 

2fc^ [\g k (n;v)\ 2 - | 5 fe(«-fc;?7)| 2 ] P n (t) 



n=0 



Defining the coefficients Fk(n; 77) by the relation 
V 2k \K\ 2 F k (n; V ) = k [\g k (n; V )\ 2 - \g k (n-k; V )\ 2 
one obtains for Eq. JT4] ) 
d\n{r)) 



dr 2 



with the (dimensionless) scaled time r given by 
t = V2r) k \K,\t. 



(14) 



(15) 



(16) 



(17) 



The coefficients F k (n]rf) determine the motional dynam- 
ics and especially the existence of a divergence in finite 
interaction times, as depicted in Eq. (^|). From Eqs. (||) 
and (|l5| ) the coefficients follow as 



Fk(n\rj) = k 



(n+k)l 2 

71 < 



(n~k) 



(18) 



with k > 1 and the functions /fc(n;7/) given by Eq. 
with fk(n;rj) = for n < 0. 



IV. LAMB— DICKE APPROXIMATION: 
EXPLODING SOLUTIONS 

In the Lamb-Dicke approximation (77 — ► 0) the coeffi- 
cients Fk(fi;rf) read as 



F k (n-0) = 



1 



(fc-1) 



n + k 
k 



(19) 



From Eq. ( |l9| ) it can be seen that, in general, the func- 
tions F k (n; 0) are polynomials in n of the order k— 1, that 
is 



fc-i 



F k (n;0) = J2 



(iki n\ 



(20) 



with a non- vanishing highest-order coefficient a k , k -i 7^0. 
From Eq. (|l9| ) it can be seen that the expansion coeffi- 
cients are always positive a k i > 0. Moreover, the lowest- 
order coefficient a k o is non-vanishing, since 



^0 = ^(0:0) 



1 



(fc-l)! 



> 0. 



(21) 



The second-order differential equations for the mean 
excitation number in the Lamb-Dicke approximation, 
Eq. (Ill) together with Eq. (|20|), read as 



d 2 (n 



dr 



( # = |> (n l (r)). 
1=0 



(22) 



Note, that the second-derivative (|22j) is always positive 
and non-zero due to the non- vanishing lowest order term 
afco, cf. Eq. (pil"|). For obtaining a lower bound for the 
second-derivative we use the relations following from the 
Schwarz inequality, 



(fi'(r)) > (n(r))', 2 = 0, 1,2,. 



(23) 



Due to the positiveness of the coefficients a k i a lower 
bound of the right-hand side of Eq. (^2|) follows by using 
Eq. P), 



d 2 (h(T)) ^-1 , 

1=0 1=0 



dr 2 



(T))>^2a k i{n(r)) l >0. (24) 



From the formal solution of Eq. ( p2] ) and the application 
of Eq. (HJ) one obtains the inequality 



(n(r)) = n + n' Q r + / dr' / dr 

Jo Jo dr 



„ d 2 (n(r")) 



(25) 



> n + rin r 



/ W ^"E« 

Jo Jo , =0 



(n(r")) ( , 



with the initial conditions 



no = («(t))| t=0 , n = 



d(ft(r)) 



dr 



(26) 



r=0 



where we have chosen, without loss of generality, the ini- 
tial time to be r = 0. From Eq. (^||) it follows, that the 
solution N\b(r) of the second-order differential equation 



d 2 N lb (r) 



k-l 



= $> w <(t). 



dr' 

L — U 

obeys the relation 

JVib(r) < (n(r)), (r > 0), 



(27) 



(28) 



i=0 



for r>0 and identical initial conditions iVn>(0) = no and 
N{ h (0) — n' . That is, M b (r) represents a lower-bound 
(lb) for the solution (n(r)). We will show in the follow- 
ing, that for k > 3 the lower-bound solution iVib(T) may 
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diverge for finite interaction times, so that it is proved, 
that the correct solution (n(r)) also diverges. 

We obtain a first-order differential equation by consid- 
ering the first derivative N[ h (T) = dN\ b (r)/dT, 



dN( h 
dr 



dN[ h dN lb 
dN lh dr 



dN' 



fc-i 



115 



d/V, 



lb 



(29) 



z=o 



The last equality in Eq. 
integration of 



J29) can then be easily solved by 



fc-i 



1=0 



N[ b dN{ h = a ki N} b dN lb , 

and one obtains 
~dN lb ( T ) ' 



dr 



(30) 



(31) 



with the coefficients bki — 2a,k,i-i/l > and bki ^ 0, 

To demonstrate the unphysical properties of the 
Hamiltonian Hj^ D ^ for k > 3 in the Lamb-Dicke approx- 
imation (|5|), it is sufficient to prove the unphysical be- 
haviour for one physically reasonable initial condition. 
For the special cases k = 3, 4 this has been already ex- 
plicitly shown in Refs. here we want to show the 
unphysical behaviour, in a general way, for all fc>3. For 
the atom initially (at r = 0) in its vibrational ground- 
state, 1^(0)) = |0), the initial conditions are no = n' = 



[the latter can be seen from Eq. ( |1C| ) with ip n (0) =5 n ,o\ 
and the differential equation j3l| ) reduces to 



dN lh (r) 
dr 



U=l 



(r) 



1/2 



(32) 



Here we have chosen the positive square-root since for 
vanishing initial velocity, n = 0, and always positive ac- 
celeration, d 2 Ni h (T)/dr 2 >0 [cf. Eqs. (§J) and @], the 
velocity at time r>0 has to be positive, dN\ b (T)/dT>0. 
Eq. j32| ) can then be integrated from the finite time t\ > 
to T2 > T\ which gives the relation 



Ti - T\ 



Mb(T 2 ) 



dn 



Wi b (n) V&fcfc n k + ... + b k2 n 2 + b kl 



(33) 



Since the velocity dN\ b {r) / dr is always positive and non- 
vanishing for t > 0, it is clear, that N\ b (Ti) > no = 0. 
Therefore the integration starts with a positive and non- 
vanishing value of the excitation, N\ b (ri) > 0, that has 
been attained after the interaction time ti . 

Now we are interested in the further evolution in the 
time interval ti—t\. In particular, we are looking for 
that time interval At^ = t 2 — t\ for which the excitation 
number N\ b (r^) attains an infinite value N\\,(t^) — > oo. 
By taking only the highest-order term in the square-root 
of Eq. we obtain an upper bound for Ar^ 



At™ < 



dn 



oo, 



ATib(Ti) 



V b kkn k 



k-2 



^b fcfe < b - 2 (Ti)' 



(fc = l,2), 
(fc>3). 

(34) 



This result reveals that we obtain a finite value of At,^ 
for k > 3. That is, after attaining the finite (non- 
vanishing) excitation N\ b (ri) after the interaction time 
Ti, the solution Ni b (r) of Eq. (27) diverges already af- 
ter the finite time interval At^, according to Eq. (p4). 
Concluding, for fc>3 the solution of Eq. (|2|) [(n(T 2 )) > 
Nihfa)] will also diverge at a certain finite interaction 
time T2 < ti + Atoo. For k — 1, 2 no upper bound for the 
interaction time can be given and it can be seen by direct 
integration of Eq. (22) that the mean excitation does not 
diverge in a finite interaction time. 



V. OVERLAP OF MATTER AND LIGHT WAVES: 
REGULAR BEHAVIOUR FOR LARGE 
EXCITATIONS 

In the preceding section it has been shown that in the 
Lamb-Dicke approximation the mean motional excita- 
tion number diverges in a finite interaction time for the 
cases fc>3. In this section we will proof, that the exact 
Hamiltonian ([!]), e.g. without the Lamb-Dicke approxi- 
mation, does not exhibit such a divergence problem. This 
is due to the overlap of matter and light waves described 
by the nonlinear operator functions (|2|). They lead to an 
excitation-dependent coupling strength which suppresses 
the unbounded increase of the mean excitation. 



A. Proof of the regular behaviour 

To prove the regular behaviour of the dynamics of the 
system described by Eqs. ([!(]) and ([l8]) we may consider 
the following situation: 

(a) If the mean quantum number would diverge, we would 
be operating in a regime of very large quantum numbers 
n. Therefore we are allowed to use an asymptotic expan- 
sion of the coefficients Fk{n;rf) for large n. 

(b) Since for k = 1, 2 we know that in the Lamb-Dicke ap- 
proximation [described by Fi^(n;0)] the dynamics does 
not exhibit a divergence in finite interaction times, it 
is sufficient to show that the asymptotic expansion of 
Fk(n; rf) has an upper bound leading to a dynamics which 
is at least as convergent as for Fx^{n] 0), 



F k (n;r)) < F lt2 (n;Q), (n » 1). 



(35) 



Then the acceleration d 2 {n(r))/dT 2 is always smaller 
than those for the well-behaved cases and a divergence in 
finite times cannot exist, regardless of the initial motional 
quantum state chosen. 
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We start by expressing the function F k (n;rj) given 
in Eq. ( fl8|) in terms of Laguerre polynomials by using 
Eq. ©, 



F k {n;rj) 



(n + k)\ 
(n-k)\ 



(36) 



V/2 



While the first (positive) term in Eq. ( p6| ) describes the 
transition to higher- lying states \n) —T\n + k), the sec- 
ond (negative) term describes transitions to lower-lying 
states \n) —>\n—k) , leading to a decrease of the accelera- 
tion. An upper bound for F k (n\ rj), which determines the 
maximum acceleration, is therefore given by neglecting 
the transitions to lower-lying states (which do not cause 
a divergent behaviour) 



F k (n;r}) < 



(n + fc)! 



4 fc V; 



V/2 



(37) 



Using the relation between the Laguerre polynomials 
and the confluent hypergeometric (Kummer's) function 

M(a,6;x)H 

L^(x)=h^M(-n,k + l;x) ) (38) 
one arrives at the inequality for F k (n; 77) 

F k (n; V )<^ n+ k k ^j M 2 (-n, k + l; V 2 ) e^ 2 . (39) 

An asymptotic expansion of the confluent hypergeomet- 
ric function M(a, 6; x) for a—t— 00, bounded b, and real- 
valued x is given by p3| , 



M(a,b;x) ~ 



1 i-ib 



(40) 



X cos 



J(2b-Aa)x - -bn + -it 
2 4 



Thus, for large numbers n the inequality reads in its 
asymptotic form 



F k (n;r]) < 



1 (n + k)\ 



> ; 1 + fc 

V n-\ — 



1 + fc 1 1 
2ri\jn+— 5 4^ 



(41) 



Therefore, the maximum value of the right-hand side of 
Eq. ( f4l| ) can be estimated by taking the squared cosines 
to be unity, 



F k (n;n) < 



(n + k)\ 



.> / 1+fc 



(42) 



Moreover, the expression ( Egfl can be further estimated 
by the following relation 



(n + fc)! , , . . 

* -L = (n + fc (n + fc - 1 . . . (n + 1 

n! 

< (n + fc) fc , 
which gives one a further simplification, 



F k {n;rj) < -e^' 2 



V 2 (n + ^f) 



(n + fc) 
» 2 (n+^) 



(43) 



(44) 



For the range of large numbers n, we are considering here, 
the function therefore has the following upper bound, 



1 e^l 2 1 
F k {n\ri) < —=, 



(45) 



that is, for large numbers n the upper bound of the func- 
tion F k (n; rf) decays as 1/y/n. It therefore can be further 
estimated by a simple constant C k {rj), 



F k (n;r,)<C k (r,), C{r,) = - 



(46) 



The resulting differential equation for the upper bound 
(ub) A^ u b(r) of mean quantum number reduces then for 
possibly large numbers n to 



d 2 N uh (r) 
dr 2 



= C k {ri). 



(47) 



Reconsidering the formal solution (|25| ) and the upper- 
bound acceleration (^) it becomes clear that A r u b(r) in- 
deed is an upper bound for the exact mean excitation 
number, 



N uh (r) > (n(r)>, 



(48) 



for identically chosen initial conditions A r u b(0) = no, 
N^ h (0) = n' o and large excitations, (n(r)) > 1. Eq. @ 
states, that the mean excitation number does not di- 
verge in finite time, since the differential equation for 
large values of n leads to a behaviour which is as con- 
vergent as in the case of fc = 1 in the Lamb-Dicke limit 
where Fi(n;0) = 1 = const. That is, the upper-bound 
solution of Eq. ([47|), N uh (r), which can be obtained by 
direct integration, 

JV ub (r) + V + (49) 

does not diverge for finite interaction times r. In conclu- 
sion, it has been proved that the mean motional excita- 
tion number resulting from the full Hamiltonian (Q) does 
not diverge for finite interaction times. 
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FIG. 2. Exact time evolution of the mean motional exci- 
tation number (n(r)) for k = 3 and Lamb-Dicke parameter 
r] = 0.2, as a function of the scaled time r given in Eq. (0). 



B. Numerical examples 

As an example, we show in Fig. || the exact time evo- 
lution of the mean motional excitation number (n(r)) 
for k — 3 and Lamb-Dicke parameter r\ = 0.2. It clearly 
shows, that instead of diverging in a finite interaction 
time as would be expected in the Lamb-Dicke approx- 
imation, the mean excitation number exhibits an oscil- 
latory behaviour. This is due to the destructive overlap 
of matter and light waves leading to a decoupling of the 
atomic motion from the laser fields for certain excitation 
amplitudes. 

To gain more insight into the distribution of the quan- 
tum state in phase space, we show in Fig. [| the time evo- 
lution of the Q function for the three-quantum coupling 
(k = 3) and for rj = 0.2. It can be seen that the dynamics is 
strongly modified by the occurrence of the circles of van- 
ishing coupling strengths. In contrast to the dynamics 
in the Lamb-Dicke approximation, where the "star" -like 
structure would be extended to infinitely large phase- 
space amplitudes, the extension of the "star" structure 
is halted at the first circle of vanishing coupling. Parts of 
the phase-space distribution are smoothed over the circle. 
For those components of the distribution that accumulate 
a phase shift of 7r/3 relative to the initial "star" struc- 
ture, the Hamiltonian (|l|) effectively exhibits a change 
of sign accompanied by a reversal of the time evolution. 
Consequently, those components of the quantum state 
are moving back towards the origin of phase space |24| . 
This effect explains the decrease of the mean motional 
excitation number as seen in Fig. |[ Note that the distri- 
bution in Fig. |3|(f) for time t = 5.74 corresponds to a local 
minimum of (n(r)) in Fig. |^. Obviously, there are some 
components of the phase-space distribution which cross 
the barrier. However, because of the existence of further 
barriers at approximately equidistant radii, the explosive 
dynamics occurring in the Lamb-Dicke approximation 
and also in the optical parametric approximation does 
not occur. 



VI. SUMMARY AND CONCLUSIONS 

In summary it has been shown that for a trapped atom 
which is driven by Raman-laser fields, in the Lamb-Dicke 
approximation a behaviour appears which is analogous to 
the case of /c-photon down conversion in nonlinear optics. 
A divergent behaviour of the mean motional excitation 
number after finite interaction times occurs for higher- 
order quantum couplings with fc>3, similar to the situ- 
ation for the parametric approximation in nonlinear op- 
tics. We have discussed these divergences within a single 
unified framework for all orders k > 3. Moreover, it has 
been argued that the Lamb-Dicke approximation, which 
is only valid for well-localised atoms, is not consistent 
with the occurrence of large (or even diverging) mean 
excitations. 

To overcome the divergent behaviour, one has to treat 
the full problem without the Lamb-Dicke approxima- 
tion. This includes the correct description of the laser- 
induced momentum transfer onto the center-of-mass of 
the trapped atom. These are described by a nonlinear 
operator function, which plays an essential role for the 
dynamics of the motional quantum state of the atom. 
By using an asymptotic expansion, we have proved that 
the correct description of the recoil effects widely modi- 
fies the dynamics for large excitations and prevents the 
mean excitation number from exploding for finite inter- 
action times. That is, the full problem leads to a regular 
dynamics where the energy of the motional degree of free- 
dom does not unphysically diverge. On the other hand, 
the Lamb-Dicke approximation fails for these types of 
couplings, as does the parametric approximation in non- 
linear optics. Whereas in nonlinear optics the divergence 
problem arises from the neglection of the pump-mode de- 
pletion and entanglement of the involved field modes, in 
the case of a trapped atom the unappropriate treatment 
of the recoil effects in the Lamb-Dicke approximation 
leads to the unphysical behaviour. 
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